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I. Introduction 



The objects of our study are the densest packings, particularly of balls and 
polyhedra, in a space of infinite volume. These optimization problems are of a 
fundamental nature and have been of interest for many years; for a survey see the 
classic texts [Feje] and [Roge], and the review [FeKu]. They appear as part of 
the eighteenth problem in Hilbert's list [Hilb], a problem devoted to fundamental 
domains in spaces of constant curvature, and the manner in which the domains fill 
up such a space. 

Most interest has centered on densest packings in the Euclidean spaces E"^, 
notably when the dimension n is 2 or 3, but we will see that analysis of the problem 
in hyperbolic spaces can clarify some issues for the Euclidean problems so we 
consider the more general problem in the n dimensional spaces where §" will 
stand for either E" or H"'. (It would be reasonable to generalize our considerations 
further, to symmetric spaces, and even to include infinite graphs, but as we have 
no noteworthy results in that generality we felt it would be misleading to couch our 
considerations in that setting.) 

Let us begin with some notation and basic features of density. Given some 
finite collection B of "bodies" in - a body being a connected compact set with 
dense interior and boundary of volume - we consider "packings" of by the 
bodies: that is, collections P of congruent copies of the bodies in which interiors 
of bodies do not intersect. Denoting by B^ip) the closed ball in of radius r and 
center p, we define the "density relative to Br{p)" of a packing P as: 

where m§" is the usual measure on Then, assuming the limit exists, we define 
the "density" of P as: 

It is not hard to construct packings P for which the limiting density D(P) does 
not exist, for instance by the adroit choice of arbitrarily large empty regions so that 
the relative density oscillates with r instead of having a limit. (In hyperbolic space 
the limit could exist but depend on p, which we also consider unacceptable.) The 
possible nonexistence of the limit of (2) is an essential feature of analyzing density in 
spaces of infinite volume; density is inherently a global quantity, and fundamentally 
requires a formula somewhat like (2) for its definition [Feje], [FeKu]. We discuss 
this further below. 

Certainly the most important examples for which we have significant informa- 
tion about the densest packings are the densest packings of balls of fixed radius in 
E" for n = 2 and 3. (For a recent survey of this problem in higher dimensions see 
[CoGS]). It will be useful in discussing these problems to make use of the notion of 
"Voronoi cell" , defined for each body /3 in a packing P as the closure of the set of 
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those points in closer to /3 than to any other body in P. A noteworthy feature of 
the n = 2 example is, then, that in the optimal packing (see Figure 1) the Voronoi 
cell of every disk (the smallest regular hexagon that could contain the disk) has the 
property that the fraction of the area of this cell taken up by the disk is strictly 
larger than for any other Voronoi cell in any packing by such disks. (Intuitively, 
the optimal configiuation is simultaneously optimal in all local regions.) As for 
n = 3, it is generally felt that the densest lattice packing (i.e. the face centered 
cubic) achieves the optimum density among all possible packings, along with all the 
other packings made by layering hexagonally packed planar configurations, such as 
the hexagonal close packed structure; see [Roge] . There are claims in the literature 
by Hsiang [Hsia] and by Hales [Hale] for proofs of this, and there is hope that the 
problem will soon be generally accepted as solved. 

Less well known but perhaps next in significance as examples of optimal density 
(see [Miln]), are the various "aperiodic tilings", especially the "Penrose kite & dart 
tilings", the tilings of by congruent copies of the two polygons of Figure 2. (A 
portion of a kite & dart tiling is shown in Figure 3.) A key feature of these bodies is 
that the only way to tile the plane with them is with a tiling whose symmetry group 
does not have a fundamental domain of finite volume; this situation is the defining 
characteristic of "aperiodicity" , and has led to renewed study of the symmetry 
of tilings (and thus packings); see [Radi]. There are other significant symmetry 
features of this example which will be discussed below, where they will help to 
develop an appropriate notion of equivalence among optimally dense packings. 

In practice it is almost impossible to actually determine optimally dense pack- 
ings - for instance, there is nothing yet proven qualitatively of the densest packings 
in the Euclidean plane by regular pentagons of fixed size (see however [KuKu]) - 
and we follow the lead of Hilbert and others in concentrating on general features 
of optima, such as their geometric symmetries. (It is because of their qualitative 
symmetry features that we attributed high significance to the aperiodic tilings.) We 
will discuss below some optimization results for packings in hyperbolic spaces, by 
balls in H", and by a certain polygon in H^. 

We note that, for optimization in a Euclidean space, there is no difficulty of 
existence of an optimum density, even though, since the limit of (2) does not exist 
for some packings, we are not able to make a comparison among all packings. One 
way to understand this is to use the fact that the relative densities dB^{p) (P) of any 
packing P can be well approximated by packings P' that have compact fundamental 
domains; intuitively, the supremum of the densities of such symmetric packings P' 
is the desired optimum density, and can be shown to be achieved, in the sense of 
(2), by some packing which is a limit of such symmetric packings. (See [BoR2] for 
a complete argument.) 

The situation for hyperbolic space packings is much more complicated, due to 
the fact that the volume of a ball of radius R grows exponentially in R ([Bear], 
[Kato]). This has the consequence that a significant fraction of the volume is near 
the surface of the ball, so it is by no means clear that one could make a useful 
approximation using a packing with cofinite symmetry group. Once one is prevented 
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from reducing the problem to such symmetric packings, one is confronted with the 
difficulty of showing the existence of a limit such as (2) for a purported optimal 
packing; see [BoR2] for a history of this difficulty. In summary, lack of proof of the 
existence of appropriate limiting densities was an impediment to progress in the 
study of optimal density (and, a fortiori, optimally dense packings) in hyperbolic 
spaces for many years [BoR2]. This led to a search for alternatives to the notion 
of optimal density. The best known of these are that of "solid" packings, and 
"completely saturated" packings [FeKu]. Both notions are defined through local 
properties of packings. The quest for a local approach/alternative to optimal density 
is perhaps reasonable given that the only practical way yet devised to prove that a 
packing is optimal is to prove that it is locally optimal in all local regions (as noted 
above for disk packing in Ep, and as the basis for the methods for ball packings in 
[Laga]). However, no local alternative has proven satisfactory [FeKu], [Bowe]. 

So far we have concentrated on the question of existence of solutions (i.e. opti- 
mally dense packings) for our general optimization problem, especially the difficul- 
ties for packings in a hyperbolic space. This existence problem was solved recently 
([BoRl], [BoR2]), where the main obstacle, the existence of limiting densities, was 
obtained through an ergodic theory formalism outlined in the next section. One 
goal of this paper is analysis of the uniqueness of solutions for our optimization 
problems. Currently, there are difficulties even for Euclidean problems. For in- 
stance, even though it would be intuitively satisfying to declare that the problem 
of optimally dense packings of by disks of fixed radius has the "unique" solution 
discussed above (Figure 1), there has been no satisfactory way to exclude some other 
packings of the same density, for instance those obtained by deleting a finite number 
of disks from this packing [CoGS]. This has been a serious obstacle to treatment 
of optimal density as one treats other optimization problems [Kupe] , and will be a 
useful guide for our approach to a general understanding of the qualitative features 
of optimally dense packings. 

The situation is significantly more complicated, and interesting, for the optimal 
packings (tilings) by kites Sz darts in (Figure 2) than it is for equal disks. It 
can be proven that there are uncountably many pairwise noncongruent such tilings, 
but that every finite region in any one tiling appears in every other such tiling, so 
they are, in some sense, "locally indistinguishable" [Gard]. It is natural to want to 
declare that this optimization problem also has a "unique" solution, and this has, 
in effect, been the practice of those studying aperiodicity, a practice we will follow. 

The notion of uniqueness must differentiate between the situations for densest 
packings of equal balls in E^ from that of E^; for the former the expected solution 
is intuitively far from unique, containing for instance the face centered cubic and 
also the hexagonal close packed structures, which must be considered different if the 
notion is used at all. However at a deeper level it must also give a useful criterion 
for when two such optimal structures are "the same", that is, it must give some 
useful notion of the geometric symmetry of optimally dense packings, for instance 
of the kite & dart tilings. In a later section we will discuss the connection between 
our approach and the use, by Connes and others, of noncommutative topology to 
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understand the symmetry of structures such as the kite & dart tihngs; see the 
expository works [Conn], [KePu] and references they contain. 

As we hope to demonstrate, study of the uniqueness problem for packings of 
hyperbohc space wiU be useful even for understanding packings in Euclidean space. 
Much of this paper consists of the analysis of a specific family of examples of optimal 
density in a hyperbolic space, the first aperiodic examples in hyperbolic space for 
which explicit optimally dense packings have been determined. (See [BoR2] for 
an analysis of other examples of tilings in H^, such as those in [MaMo], [Moze] 
and [Good]). These examples exhibit features not seen in Euclidean examples, and 
will help us draw some conclusions about the general features of optimum density 
problems. 

II. The Ergodic Theory Formalism for Optimal Packing 

One concern of this paper is with the uniqueness of solutions to optimal den- 
sity packing problems. While proving existence of solutions for problems set in 
Euclidean spaces did not require giving the problems a formal structure, the ques- 
tion of existence of solutions for problems set in hyperbolic space definitely did 
require introducing a formal structure, and we will see that this same structure is 
useful for handling questions of uniqueness, even in Euclidean space. We follow 
[BoRl], [BoR2] in introducing an ergodic theory structure into our optimization 
problems, in order to control the existence of limits such as (2). 

Using the notation of section I, consider the space Vb of all possible packings 
of by bodies from B, and put a metric on Vb such that convergence of a sequence 
of packings corresponds to uniform convergence on compact subsets of S"^. Such a 
metric makes Vb compact, and makes continuous the natural action on Vb of the 
(connected) group Q'^ of rigid motions of [RaWo] . 

We next consider Borel probability measures on Vb which are invariant under 
G^- (To construct such a measure consider any packing P for which the symmetry 
group has fundamental domain of finite volume, and identify the orbit 0{P) of P 
under as the quotient of Q"' by that symmetry group. One can then project Haar 
measure from Q'^ to an invariant probability measure on 0{P) and then extend it 
to all of Xb so that the complement of 0{P) has measure zero.) 

We define the "density of the (invariant) measure /x" on Vb, D{fi), by D{fi) = 
IJ,{A), where A is the following set of packings: 

A = {PeVB\ the origin O of is in a body in P}. (3) 

(It is easy to see from the invariance of fi that fi{A) is independent of the choice of 
origin.) We may now introduce the notion of optimal density. 

Definition 1. A probability measure on the space Vb of packings, ergodic under 
rigid motions, is "optimally dense" if D{fi) = sup^ D{fx) = sup^yu(A); the number 
sup^ H{A) is the "optimal density" for packing bodies from B. 

(An invariant measure n is "ergodic" if it cannot be expressed as an average: 
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ji = aiHi + a2//2) with oi, 02 > and jii ^ 112 invariant.) It is not hard to show 
[BoRl] the existence of such optimal measures for any given B. 

Finally, the terminology is justified as follows. First we rewrite the right hand 
side of (1) as: 

where Xa is the indicator function for A, v is Haar measure on Q'^ and 

g^{r,p) ^{geg^\ dsr.[g{p),p\ < r}, (5) 

where dsn is the distance function on It follows from G.D. Birkhoff 's pointwise 
ergodic theorem [Walt] that for = E"' and any crgodic fi there is a set of P's, of 
full /U-measure and invariant under Q'^, for which the limit in (4) exists. This has 
been extended to S"^ = by Nevo et al: [Nevo] , [NeSt] (with the invariance of the 
set of P's proven in [BoR2]). We may conclude then that "most" of the packings 
in the support of a fixed ergodic measure have the same well defined density in the 
sense of (2) ; so as one varies the measure one sees Vb decomposed into packings 
of various densities, with those of optimal density being the ones in which we are 
interested. Formally we define optimally dense packings (slightly more stringently 
than in [BoR2]) as follows. 

Definition 2. A packing P is "optimally dense" if it is "generic" for some optimally 
dense /i, that is, it is in the support of fi and: 
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fiQ) dM) = lim ^^^^r-^\ [ f^siP)] di^ig), (6) 



for every 7? e and every continuous / on Vb- The set of all optimally dense 
packings for bodies in B will be denoted Pg. 

We note that the set of packings generic for the invariant measure fj, is of full 
measure with respect to this follows from the ergodic theorem and the fact that 
the space of continuous functions on Vb is separable in the uniform norm. 

Lemma 1. If P is generic for the ergodic measure /i (which is not necessarily 
optimally dense) then: 

lim rrn) M / XA[g{P)]dv{g)= ! XA{Q)di,{Q) 

r^oo u[g''{r,p)] J gr.(^r,p) Jvs (7) 

= D{^i) 

for every p G 
Proof. Let 

A' = (P e Vb | C is in the interior of a body in P}. (8) 
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Define continuous fk on Vb by: 



{1, on A 
0, on{PeA-\dsr.{0,dP)>^} (9) 
1-kc, on{PeA'=\d§n{0,dP) = c<l}, 

where dP denotes the union of the boundaries of the bodies in P. Note that the fk 
decrease pointwise to Xa- Similarly define continuous gk on Vb by: 

rO, on A'" 

gk{P)={h on{PeA'\ds40,dP)>^} (10) 
[ kc, on {P e ^' I dsr. (O, aP) = c < i}. 

Note that the gk increase pointwise to Xa'- Now, given e > 0, choose K > such 
that 



< J fKdfJ,- J XAdfj, < e/2 and < J Xa' dfj, - J gxdfj, < e/2. (11) 
Define, for i? > 0, the measure i^(i2, P) on Vb by 

for continuous /, where as before u is Haar measure on Q"^. Then choose R > 
such that 

I J fKdu{R,P) - J fKdfi\ < e/2 and | J gKdu{R,P) - J gKdfi\ < e/2 (13) 



for all R> R. We then have: 



JxA'dii-e< J gKdiy{R,P) < J Xa' diy{R, P) < JxAdiy{R,P) 



< 



J fKdv{R,P) < JxAdfi + 



(14) 



However, from the ergodic theorem J Xa—Xa' dfi = J Xa/a' dfi = 0, so J Xa' dv[R, P) 
and J Xa du{R, P) both converge to J Xa dfi = D{p) as R oo. | 

Next we note a useful tool for computing optimal densities. For those P G Vs 
such that the point p & is contained in the interior of a Voronoi cell (which cell we 
denote by Vp{P)), we define Fp{P) to be the relative volume of Vp{P) occupied by 
the bodies of P. (We note that Fp is defined //-almost everywhere for any invariant 

Definition 3. For invariant measures /i we define the "average Voronoi density for 
Dvipt), as jj,^ Fp{P) dii{P). (Note that Dv{lJi) does not depend on p because 
of the invariance of ji.) 
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The notion of average Voronoi density is useful, as it has been shown [BoRl, 
BoR2] that, for any invariant measure ji, the average Voronoi density Dy^ij) equals 
the average density D^jj). 

Summarizing the above, we have sketched a formalism through which one 
proves existence of solutions to the general problem of densest packings of by 
congruent copies of bodies from B. Our next goal is to consider uniqueness, but 
we will first need to discuss symmetry further (section III), and a new family of 
examples (section IV). 

III. Symmetry in the Problem of Optimally Dense Packings 

We will use the following common terms: a packing is called "periodic" (rcsp. 
"nonperiodic" ) if it has (resp. does not have) a symmetry group with fundamental 
domain of finite volume, and we say an optimal packing problem is "aperiodic" if 
all its optimally dense packings are nonperiodic. 

Although it will be more difficult to introduce a general criterion for uniqueness 
(which we will attempt in section V), the above formalism can be used to solve the 
old problem of making sense of uniqueness for the problem for disks of fixed radius 
in E^. 

Theorem 1. There is only one optimally dense packing in E^ for disks of fixed 
radius, up to rigid motion. 

Proof. Assume /x is any optimally dense measure for this problem. Using the fact 
that the hexagonal Voronoi cells of P° (the hexagonal packing of Figure 1) are, up 
to rigid motion, the unique cells of optimal density (see [Feje]), and using the basic 
result on average Voronoi cells [BoRl], we see that for |U-almost every packing P, 
the cell containing a particular point p must be this regular hexagon. Repeating this 
argument for a countable dense set of p's, we see that for //-almost every packing P 
every Voronoi cell is, up to rigid motion, this regular hexagon, i.e., //-almost every 
packing P is P°, up to rigid motion. It is shown in [BoRl] that there is a unique 
invariant measure with support in the (closed) orbit of any periodic packing, and 
therefore the orbit of P° consists of all the optimal packings. | 

The ergodic theory formalism automatically makes sense of the uniqueness of 
this optimization problem, by its subjugation of packings to invariant measures on 
packings. Aperiodic problems, such as the kites & darts, are more subtle. Our next 
result is a connection between aperiodicity and the uniqueness of packing problems. 

Theorem 2. If there is only one optimally dense packing of §"■, up to congruence, 
by congruent copies of bodies from some fixed, finite collection then that packing 
must have a symmetry group with compact fundamental domain. 

Proof. By assumption there exists a probability measure /i, invariant under Q'^, 
for which the orbit 0{P) of some packing P has measure one. Since 0{P) can be 
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identified with the quotient of by the symmetry group Tp of P, it follows that 
Tp is cofinite. Now if Tp were not cocompact (this complication is only relevant for 
hyperbolic S") it would follow (by considering cusps in S'^/Fp) that P contained 
no bodies in certain regions containing arbitrarily large balls. But it is easy to see 
that optimally dense packings must be "saturated" , that is, one cannot add another 
body to the uncovered regions of an optimally dense packing; in fact it was proven 
in [Bowe] that the set of completely saturated packings have full measure with 
respect to any optimally dense measure. (A "completely saturated packing" is one 
in which it is impossible to remove any finite number of bodies from an optimally 
dense packing and replace them with bodies of larger volume.) | 

One thing we can conclude from these two theorems is that the problems of op- 
timal density decompose naturally into two classes: those allowing periodic optima 
(such as with balls of fixed size in or E^), and the class of aperiodic problems. 
The former no longer pose any difficulty as to classifying their uniqueness, leaving 
us now to understand the more interesting class of aperiodic problems. 

Aperiodicity is not an unnatural circumstance - it may even be generic in some 
sense; see [MiRa] for a related problem. In Euclidean space aperiodicity has only 
been discovered so far in packings of complicated polyhedra, whereas in hyperbolic 
space it already appears in packings of balls. 

Theorem 3 [BoRl]. For all but countably many fixed radii R the ball packing 
problem in has only aperiodic solutions. 

For packings/tilings in Euclidean space there has not yet been a significant 
attempt to understand the uniqueness problem. We will postpone our attempt 
at a formal definition until after considering the following new examples, the first 
aperiodic problem in hyperbolic space for which we can determine explicit solutions. 

IV. A New Example: the Modified Binary in 

Using the upper half plane model of H^, consider the "binary tile" r of Figure 
4 (introduced by Roger Penrose in 1978 [Penr]). We define as the "core" 7 of r that 
shape with four edges - two segments of geodesies, and two segments of horocycles 
(one twice the length of the other) - obtained by omitting the bumps and dents of r. 
More specifically we could take the coordinates of the vertices of 7 to be i, i + 2, 2i, 
and 2i + 2. Congruent copies of r can only tile as in Figure 5. In fact it is 
useful to classify these tilings as follows. Once the location of one tile is known, the 
bumps on the geodesic edges of (the core of) r force the positions of tiles filling out 
the region between two concentric horocycles (the ones containing the horocyclic 
edges of the (core of the) tile). Consider now the possible tiles abutting the ones 
in this "horocyclic strip" . There is only one way to fill an abutting strip which is 
"further" from the common point at infinity of the horocycles, and two ways to fill 
the strip which is "closer" . This fully classifies the possible tilings of by r. 

We now construct a new tile f based on the same core 7, which will permit some 
new tilings. On the geodesic edges of 7 we add the same bumps and dents as before. 
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but we enlarge each of the other original bumps and dents as follows. Consider 
a (densest [Boro]) packing of the hyperbolic plane by horoballs, as illustrated in 
Figure 6, and consider three abutting horoballs. Divide the (white) region between 
the horoballs into 3 congruent regions by means of three geodesies, each drawn 
from the center of mass of the region to the point where a pair of horoballs touch. 
The new tile f is illustrated in Figure 7. From the construction we see that, as 
with T, in any tiling of the plane by f, once we know the location of a specific 
tile we can uniquely fill in a horocyclic strip, and then have two choices for filling, 
consecutively, each of those strips which are closer to their point at infinity, thus 
filling in a horoball. However, there is now a second way to fill the abutting strip 
which is further from the point at infinity, in which the bumps from three tiles abut 
to fill in a white region in Figure 6. If this latter method is used, the only way 
to complete a tiling of the plane is to fill in each of the horoballs defined by the 
two new tiles, then use the same method to extend beyond these horoballs to more 
horoballs, etc. Intuitively, these new tilings are obtained from a densest horoball 
packing by tiling each horoball with copies of the tile f. 

Let T(f) be the set of all tilings of the plane by f. We call such a tiling 
"degenerate" if the cores of the tiles themselves tile the hyperbolic plane, and "non- 
degenerate" if the union of the cores corresponds to a (densest) packing P of by 
horoballs. (Such a horoball packing has symmetry group conjugate to PSL(2,Z), 
which is cofinite but not cocompact.) We know from [BoRl] that the set T^eg of 
degenerate tilings has measTire with respect to any invariant probability measure 
on T(f). As a result, degenerate tilings do not qualify as optimally dense packings 
of f in the sense of Definition 2. However, invariant measures on T(f) do exist. 

To construct such a measure we consider the internal structures of the diff'erent 
horoballs. This structure is related to a choice of the two lower bumps on the tile 
f, henceforth called "prongs" . For each horoball H and each triangle that touches 
the horoball, the internal structure of H is associated to a dyadic integer, that is a 
formal sum X^^o'^i^*' with Oj e {0, 1}, where two dyadic integers are considered 
close if their first N terms agree, with N large. The first digit tells whether the left 
or right prong of a tile from H sticks into the triangle, the next digit tells whether 
that tile emerges from the left or right prong of its "parent" , and so on. We let each 
digit be an independent random variable, with equal probability of being or 1. Let 
the internal structures of distinct horoballs be independent, and be independent of 
the location of the horoballs, which is given by Lebesgue measure on iy^/PSL(2, Z). 

From the existence of an invariant measure, it follows [BoRl] that there are 
tilings in T(f) which are optimally dense packings of f in the sense of Definition 
2. It is not difficult to show from this that all optimally dense packings of f are 
tilings [BoRl]. We now see that the optimal density problem for the modified 
binary f is aperiodic: the presence of the horoballs immediately implies that the 
symmetry group of an optimal packing/tiling is at most cofinite, not cocompact, 
and, as argued in the proof of Theorem 2, a tiling by compact bodies cannot have 
a symmetry group which is cofinite but not cocompact. 

The presence of the closed invariant set Tdeg is a new feature in optimal density 
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problems. Since it is a subset of the orbit closure of every tiling, it is in the support 
of every invariant measure on T(f), even the ergodic ones. However, there is no 
invariant measure on Tdeg- This is not possible for problems set in a Euclidean 
space, since if the Euclidean group acts on a compact metric space an elementary 
fixed point argument ([Radi]) guarantees the existence of an invariant probability 
measure on that set. We will discuss this feature of T^eg further below, when we 
consider various types of conjugacy for the dynamical systems in which we are 
couching our optimization problems. 

We now turn to the construction of uniquely ergodic invariant subsets of T(f ), 
and the measures they support. Let be a dyadic integer. For each w, let 7^ be 
the closure of the class of tilings for which the sum of the three dyadics at each 
triangle is w. 

Theorem 4. 7^ is uniquely ergodic under the action of Q"^ =PSL(2, R). 

Proof. For T e %o — T^eg, there is naturally associated to T a horoball packing 
h{T). For N > Q and T G 7^ — T^eg, we let Qj^{T) be the packing obtained from 
T by removing all but the horocyclic rows of tiles closest to the boundary of any 
horoball in h{T). If T G T^eg, we let Qn{T) be the empty packing, 0. Note that 
Oat defines a continuous map from 7^ onto a compact space Vn of packings and 
that commutes with the action of Q"^. 

From Lemma 2 below it follows that Vn admits only two ergodic measures; 
one concentrated on and the other being derived from Haar measure on the 
space /H where H is the symmetry group of Qn(T) for any T & — %ieg 
([BoR2]). Since T^eg has /x-measure zero with respect to any invariant measure /x 
on Tyj ([BoR2]), the empty packing is measure zero with respect to the pushforward 
of fi. Therefore there is only one possibility for the pushforward of n, or cquiv- 
alently, the space 7^ is uniquely ergodic with respect to the cr-algebra Q^^{T,n) 
(where S^v is the Borel cr-algebra of Vn)- Since this is true for all and the 
cr-algebras ©^^(S^v) are increasing, this implies that 7u, is uniquely ergodic with 
respect to the cr-algebra (J^^ ©^^"'^(Ejv)- But 7^eg has measure zero with respect to 
any invariant measure /i on 7^ ([BoR2]) and the topology of 7^ — T^eg is gener- 
ated by Uat ©^"'^(top^) (where top^ is the topology on Vn)- So the |U-closure of 
Uiv ^jv^l^-'v) contains the Borel cr-algebra of 7^,. Thus 7^ is uniquely ergodic. It 
remains only to prove the following lemma. 

Lemma 2. For any T, T' G 7^, — Tdeg, ©at (7") has cofinite symmetry group, ©^^(7") 
is in the orbit of Qn{T') and Vn is equal to this orbit union the empty packing. 

Proof. The symmetry group of the horoball packing h{T) for T G 7^, — T^eg is 
conjugate to PSL(2,Z). The packing h{T) naturally corresponds to an infinite 
trivalent tree with the triangles of the packing corresponding to the vertices of the 
tree. We can navigate around the tree with two fundamental operations. 

If a "state" is a vertex together with a choice of one of the three edges leading 
out from that vertex, then the two operations on states are 



10 



L = Go forwards to the next vertex and bear left. 



C = Rotate countercfockwise by 120 degrees 



(15) 
(16) 



C and L obviously generate the entire symmetry group of the tree. In terms of 
PSL(2,Z), C is the elliptic element 



OT z ^ z + 1, and together they generate all of PSL(2, Z). 

Now we consider filling the horoballs with tiles f, so that to each triangle we 
can associate three dyadic integers, one for each of the horoballs that meet at the 
triangle. Different triangles that touch the same horoball will not have the same 
dyadic integer; rather, moving along the edge of the horoball counterclockwise (as 
seen from inside the horoball) will increase the dyadic number by one each step. 

If we apply the condition that the three numbers at each triangle must add 
up to w (a fixed dyadic integer), then two dyadic numbers at a triangle determine 
all the rest. If you know that a given vertex has numbers a and b then the third 
number must he c — w — a — b. One neighboring vertex has numbers b — 1 and 
a + 1, so its third vertex must be c. Each vertex determines its neighbors, and so 
determines the entire tree. Thus a tiling in standard position (that is, with a choice 
of preferred vertex and preferred edge directed out from the preferred vertex) can 
be associated to a pair (a, b) of dyadic numbers, where a is the index of the horoball 
to the left of the preferred outgoing edge, and b is the index of the horoball to the 
right. 

The action of L and C is easy to compute, namely: 




(17) 



or z 



1/(1 — 2;), and L is the parabolic element 




(18) 



L : (a, 6) — > (a + 1, c) 



(19) 



C : {a,b) (c,a), 
where c = w — a — b. It is not hard to check the following elements: 



(20) 



: (a, 6) ^ (a + 2,6- 1) 



(21) 



CL^C^ : (a, 6) ^ (a - 1,6+2) 
L^CL^C^ : (a, h)^{a + 3, 6) 
L^CL^C^ : (a, 6) ^ (a, 6 + 3). 



(22) 

(23) 
(24) 
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Now consider the effect of PSL(2, Z) on pairs (a, b) as above but taken modulo 
2-^. Since 3 and 2-^ are relatively prime, some power of L^CL^C^ sends (a, 6) 
to (a + 1,6) (mod 2^), and some power of L^CL^C^ sends (a, 6) to (a, 6 + 1) 
(mod 2-^). Thus PSL(2, Z) acts transitively on the space of pairs (a, b) (mod 2-^). 

Therefore, for aU T,T' e%u- Tdeg, Qn{T) is congruent to eAr(T'). Also, the 
subgroup that preserves the pair (a, b) (mod 2^) (and hence the first rows of 
each horoball in the tiling), is an index 2^^ subgroup of PSL(2, Z), and hence is a 
cofinite subgroup of PSL(2, M). Finally, Pn is the union of this orbit and the image 
of Tdeg, which is the empty packing. | 

We next give another property of these "fixed-sum" classes of tilings, in terms of 
dynamical conjugacy. For convenience we recall some common terms. A topological 
group G acts continuously on a compact metric space X if there is a map </> : 
{g,x) ^ G X X —>■ g{x) e X which is continuous and satisfies h[g{x)] — [hg]{x) for 
all g,h G G and x E X. Assuming G acts continuously on X and Y, the actions 
are called "topologically conjugate" if there is a homcomorphism a : x E X ^ 
a[x] e F such that a [(/(a;)] = (/(ctfa;]). Assume further the existence on X and Y of 
Borel probability measures fix and fiy which are invariant under the corresponding 
actions of G. These two actions are said "mcasiuably (or metrically) conjugate" 
if there are invariant subsets Xq C X and Yq C Y, each of measure zero, and 
an invertible map a' : x E X/Xq a'[x] G Y/Yq such that q;'[(7(x)] = (7(a'[x]) 
which, together with its inverse, is measure preserving. Finally we introduce an 
intermediate form of conjugacy (related to "almost topological conjugacy" [AdMa]) 
as follows. The actions of G on X and Y will be called "almost conjugate" if there 
are invariant subsets Xq C X and Yq G Y, each of measure zero with respect to 
all invariant measures, and a homcomorphism a' : x E X/Xq — > a'[x] G Y/Yq such 
that Q!'[5f(a;)] = g{a'[x]). 

Theorem 5. 7^ and 7^' are topologically conjugate if and only if w — w' G 3Z. 
7^ and 7^/ are almost conjugate for any w, w' . 

Proof. If — It; = 3e, where e G Z, then we construct a conjugacy by leaving 
the location of all the horoballs fixed, and simply adding e to the dyadic index of 
each horoball. In the A'"-th layer of a horoball the conjugacy essentially acts by 
translation by e/2^, so points deep within a horoball are moved only slightly. In 
the case of degenerate tilings, the conjugacy leaves the entire tiling fixed. 

li w — w' is not 3 times a rational integer, then {w — ?«')/3 is still a dyadic 
integer, since 3 is a unit in the ring of dyadic integers. Adding {w — w')/2t to the 
index of each horoball is a continuous map on the complement of Tdeg, but is not 
uniformly continuous and does not extend to all of 7^. This shows that 7^ and 7^/ 
are almost conjugate for any w,w'. 

The proof that w — w' G 3Z is necessary for topological conjugacy is harder, 
and will consist of four lemmas. 

Lemma 3. Any topological conjugacy between 7^ and 7^/ must preserve the 
points on the sphere at infinity that are tangent to horoballs. Furthermore, the 
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"radii" of the horoballs can only change by a finite amount. That is, there exists a 
constant R (depending only on (p) such that, if T G 7^ is a tiling and if is a horoball 
in h(T), and H' is the corresponding horoball in h[(p{T)] (that is, with the same 
tangent point), then H is contained in an i?-neighborhood of H' and vice- versa. 

Lemma 4. The R of the previous lemma is actually zero; topological conjugacies 
preserve the locations of horoballs exactly. 

Lemma 5. Let H be any horoball in h{T) for any T E T^, and let a be its 
dyadic index (measured from a particular triangle). Let a' be the dyadic index of 
H' C /i[(/)(T)] measured from the same triangle. Then the set of differences a' — a 
(for all such horoballs H in all such tilings T) is a bounded subset of Z. 

Lemma 6. There is a triangle in T, with indices a, b and c, such that a' — a = 
b' -b = c' -c = {w'-w)/3. 

Proof of Lemma 3. A topological conjugacy is uniformly continuous, so if : 7^ — > 
7^/ is a topological conjugacy, then for every r' > and e > there is a radius r 
such that, if the neighborhoods of two points agree to radius r, then from uniform 
continuity and conjugacy their images agree out to radius r', up to an "e wiggle". 
More precisely, for any two tilings Ti, T2 G 7^ and points pi, P2 G H^, if there is an 
isometry of that sends a ball of radius r of pi in Ti exactly onto a ball of radius 
r of p2 in T2, then the same isometry sends a ball of radius r' of pi in (/>(Ti) to an 
e-small distortion of a ball of radius r of p2 in (p{T2), where an "e-small distortion" 
means an isometry that moves each point in the neighborhood a distance e or less. 

Now take r' to be greater than the diameter of a triangle and e to be much less 
than the diameter of a triangle. Take any tiling T for which 0(T) has a triangle 
centered at p2- We claim that T contains a triangle centered at a point pi at 
distance at most r + 1 from p2. For if not, then the r-neighborhood of p2 (call it U) 
lies completely within a horoball of h{T). But then there is a constant, f, say such 
that every ball of size r contains an r-ball such that T restricted to that r-ball is 
isometrically conjugate to T restricted to U. (Here T is thought of as the function 
from the plane to the tile f that is induced by the tiling T). This implies that 
every ball of size r + e contains a triangle of 4>{T). But this contradicts the fact 
that there are points, deep within a horoball of h{(f){T)) that are at least a distance 
f + e away from any triangle in the complement of h{(f){T)). But this contradicts 
the fact that </>(T) is made up of horoballs, some points of which are arbitrarily far 
from triangles. 

Now let H be any horoball in h{T). Since all points in H that are farther than 
r + 1 from the boundary of H are mapped into a horoball in h[4){T)] (i.e., not into a 
triangle), and since this set of points is connected, H lies within an r-neighborhood 
of a specific horoball H' C /i[(/)(T)]. This implies that H and H' have the same 
tangent point on the sphere at infinity. 

To obtain the fixed bound R, just repeat the argument for (f)~^ and take R to 
be the larger of the two constants r + 1. | 
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Proof of Lemma 4. Let R be as before. We know that preserves the location of 
each horoball, and changes its radius by at most R. The question is which horoballs 
grow and which shrink, and by how much. Consider a triangle in T, with center 
point p, where horoballs Hi, H2 and H3 meet. It is impossible for two of these 
horoballs to grow, or one to grow while a second does not change, lest they overlap. 
If two stay fixed, then the entire pattern is fixed. Thus, if there are any changes 
anywhere, then at each triangle either one horoball grows (or stays fixed) and the 
other two shrink, or all three shrink and one or more other horoballs H4, Hrj, . . . grow 
to fill up the space. There are only a finite number of horoballs within distance R of 
the triangle, so only a finite number of possible directions where the tile containing 
p in (f){T) can be pointing. By continuity there is a number such that knowing 
all tiles within the A^th collar of the triangle determines which of Hi, H2 and H^ 
grow and which shrink. In particular, knowing the first digits of the dyadic labels 
for Hi^2,3 determines which grow and shrink. In essence, all our labels should be 
counted mod 2^. 

There are 2^^ possible triples (a, b, c) of numbers (mod 2-^) that add up to 
w (mod 2^). These correspond to labels for horoballs that meet at a triangle, 
counting clockwise. For each one, either all triangles with this label have the "a" 
tile shrink, or none do. Let Sa be the set of labels for which the "a" tile shrinks, let 
Sf, be the set for which the "6" tile shrinks, and let Sc be the set for which the "c" 
tile shrinks. We will show all three sets are the whole set of triples, so all horoballs 
shrink, which is a contradiction. 

Each horoball actually meets an infinite number of triangles. By comparing 
adjacent triangles with the same horoball we see that (a, b, c) G Sa if and only if 
(a + 1, c — 1, 6) e Sa- Continuing this process, we get that either the entire orbit 
{(a + 2n, b — n,c — n)} U {{a + 2m + l,c — m — l,b — m)}, is in Sa or the entire 
orbit is out. Note that 3 is a unit in Z2JV, so we can take n = (6 + 1 — a)/3 and 
m — (6 — a — 2)/3. This means that both 



are in the same orbit as (a, b, c). If (a, b, c) ^ Sa, then in any triangle with indices 
(a + 26 + 2)/3, {a + 2b — l)/3, and (c + a — 6 — l)/3, in clockwise cyclic order, the 
horoballs with indices (a + 26 + 2)/3 and (a + 26 — l)/3 must both grow (or stay 
the same size), which is impossible taking into account the first paragraph. Thus 
(o, 6, c) G Sa- But the triple (a, 6, c) was arbitrary, so every triple is in Sa, and 
likewise in S^ and Sc- I 

Proof of Lemma 5. First we show that for each horoball a' — a must be an integer. 
For each dyadic integer x let 7!'m{x) be the fractional part of 2~"^x, and let o-m{x) 
be the integer part. The effect of adding a' — a on the mth layer of the horoball is 




(25) 



and 



a + 26-1 c + a- 6-1 a + 26 + 2 
3 ' 3 ' 3 



(26) 
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to translate the locations of the tiles by 7r^(a' — a) and to change the pattern of 
"ancestor" tiles by a^{a' — a). If a' — a is not an integer, then ^^{a' — a) docs not 
converge, so different layers deep in the horoball get shifted by different amounts, 
which contradicts uniform continuity, insofar as each piece of each layer looks like 
a piece of every other layer. 

Thus for each horoball, the difference a' — a is an integer. If these differences 
are not bounded, we can pick a sequence of horoballs Hm (with indices a^, and 
possibly in different tilings) such that 7rm.(a^ ~ ^rn) does not converge. Since with 
radius mlog(2) every patch centered on the mth layer of Hm is replicated in the 
Mth layer of Hm-, for every M > m, this lack of convergence of TTmia'm — ctm) 
contradicts uniform continuity. | 

Proof of Lemma 6. Since the differences (a' — a) take values in a finite set, the 
values of {a' — a) can be determined by knowing the first digits of (a, b, c). But 
at some triangles, it happens that a = b = c (mod 2-^) (since all allowable triples 
mod 2^ do occur, and since 3 is a unit when working mod 2^). At such triangles, 
we must have a' — a — b' — b — c' — c — {w' — w)/3 by symmetry. By Lemma 5 this 
common difference must be a (rational) integer, so w' — w & 3Z, which completes 
the proof of this lemma, and the theorem. | 

A consequence of Lemma 6, together with the fact that two adjacent horoballs 
determine the entire tiling, is: 

Proposition 1. There is a unique topological conjugacy from 7^ to 7^/ when 
w — w' & 3Z. Equivalently, there are no nontrivial automorphisms of 7^. 

Although the tilings in T(f) - in particular those in any %u - cannot have a 
cofinite symmetry group, those in any 7^ in fact do have a nontrivial symmetry 
group, as we see next. 

Let To e 7^ — Tdeg and consider the set h~^[h{To)] of nondegenerate tilings 
having the same associated horoball packing as Tq. We shall describe the symmetry 
group of h~^[h{To)], i.e., the subgroup of PSL(2,R) consisting of those elements 
which fix every tiling in h~^[h{TQ)]. Up to conjugacy, this group is independent of 
w and Tq. 

The symmetry group of the horoball packing h{To) is conjugate in PSL(2, M) 
to PSL(2, Z); choosing a particular conjugacy is the same as choosing a triangle and 
a distinguished vertex in /i(To), i.e., a state in the trivalent tree. Fix a conjugacy 
and identify the symmetries of h~^[h{To)] with the elements of PSL(2, Z) which fix 
every pair of dyadic numbers. We recall some terminology from Lemma 2 concerning 
particular elements of PSL(2, Z). 

Theorem 6. The elements and = C^LC^L freely generate a subgroup £ 
of PSL(2, Z) of index 6. The symmetry group of h~^[h{TQ)] is the kernel of the 
abehanization (L^, i?^) i-^- Z © Z. 
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Proof. The actions of and on pairs of dyadics are given by 

i?2 : (a,6) ^ (0 + 1,6-2) and : (o, 6) ^ (a + 2, 6 - 1). (27) 

One readily checks that the operations , , , , and L^R~^ are 

precisely the ones which move from a vertex in the trivalent tree to a vertex two 
edges away and induce maps of the form (a, 6) 1-^ {a + k,b + £), k,£ & Z on pairs of 
dyadic numbers. It follows that the index of £ in PSL(2, Z) is 6. 

Freeness follows from the fact that distance from the starting point does not 
decrease as we follow some sequence of the basic operations L^^, R^^ unless one of 
the operations is followed immediately by its inverse. 

Since the vectors (1, —2) and (2, —1) are linearly independent, the symmetry 
group of h~^[h{TQ)] consists of those words in 1/=*=^, for which the sums of the 
powers of R and L are both zero, i.e., the kernel of the abelianization of £. | 

We now note that our use of the densest packing by horoballs. Figure 6 (or, 
using the Poincare disk. Figure 8), was not critical in the above method. An infinite 
family of generalizations can be made from other such horoball packings, as we now 
argue. 

To generalize our "triangular" tilings, wc consider tilings of constructed 
as follows. First pack by horoballs such that five horoballs meet along regular 
"pentagons" (rather than triangles), as in Figure 9. The symmetry group of such 
a packing, the Hecke group G5, is a cofinite subgroup of PSL(2,R) generated by 
2; — > —1/z and z ^ z + X, where A = (1 + \/5)/2 is the golden mean. 

We tile each horoball with (differently) modified binary tiles, where now we 
need an appropriate width so we can arrange that the prongs sticking out of such 
a tile each fill up a fifth of a pentagon; see Figure 10. Relative to such a pentagon, 
the tiling of a horoball with modified binary tiles is associated to a dyadic integer. 
The first digit tells whether we are on the left or right prong of the tile, the next 
digit tells whether that tile emerges from the left or right prong of its parent, and 
so on. 

Pick a rational integer k, once and for all. Let the dyadic integers around a 
pentagon, counting clockwise, be a, b, c, d and e. The "a" and "6" horoballs also 
meet at another pentagon, and wc assume that the five dyadic integers representing 
these horoballs, counting counterclockwise, are a + 1, b — 1, c — k, d and e + k. This 
rule for relating patterns around adjacent pentagons is a generalization of the "fixed 
sum" rule for triangular tilings. Notice also that the sum around the pentagons is 
fixed, and we let Tk^w be the closure of the class of tilings that follow the "fc-rule" 
and for which the sum of the five dyadics at each pentagon add up to w. Again we 
denote by Tdeg the degenerate tilings. 

Theorem 7. 7^^^, is uniquely ergodic under the action of =PSL(2, M). 

Proof. This proof is nearly identical to the proof of Theorem 4. As before, we 
approximate tilings in Tk^w by packings in which each horoball has only N layers, 
and show that such packings have cofinite symmetry groups. 
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As before let Qn{T) be the packing obtained from T by removing all but the 
horocyclic rows of tiles closest to the boundary of any horoball in h{T), and let 
Pjv be the range of Qn(T). In place of Lemma 2, we need to prove: 

Lemma 7. For any TT' G Tk^w—Tdeg, Qn{T) has cofinite symmetry group, Qn{T) 
is in the orbit of Qn{T') and Vn is equal to this orbit union the empty packing. 

Proof of Lemma 7. A "pentagonal" horoball packing of corresponds to an infinite 
5-valent tree, with the pentagons of the packing corresponding to the vertices of 
the tree. We can navigate around the tree with two fundamental operations. If a 
"state" is a vertex together with a choice of one of the five edges leading out from 
that vertex, then the two operations are 



P — Rotate counterclockwise by 72 degrees (28) 

L = Go forwards to the next vertex and bear hard left. (29) 

Together these generate G5. In terms of PSL(2,M), P is the elliptic element z — > 
1/{X — z) while L is the parabolic element 2; — > 2; + A. 

We list the horoballs around a vertex in counterclockwise order, starting with 
the one to the right of the chosen edge. We need only list the first four of the five 
horoballs, since if their dyadic integers are a, b, c, and d, then the last one must be 
e — w — a — b — c — d. 

The actions of L and P are easy to compute, namely: 

L : (a, b, c,d)^ {a + l,e + k, d,c — k) (30) 

P : (a, 6, c, d) — > (e, a, b, c). (31) 
It is not hard to check the following elements: 

: (a, 6, c, d) ^ (a + 2, 6 - 1 + fc, c - fc, d - k) (32) 

PL^P^ : {a,b,c,d) ^ {a + k-l,b + 2,c-l + k,d-k) (33) 

P^L'^P^ : (a,6,c,d) ^ {a - k,b + k - 1, c+ 2,d + k - 1) (34) 

P^L^P'^ : {a,b,c,d) ^ {a - k,b- k,c + k-l,d + 2). (35) 

Thus the possible values of (a, 6, c, d) differ by (among others) the elements of the 
sub-lattice of Z*^ generated by (2, k — 1, —k, —k), {k — 1,2, k — 1, —k), 
{-k, k-l,2,k-l) and {-k, -k, k-1,2). 
Since 



det 



^2 k-1 -k -k ^ 

k-1 2 k-1 -k 

-k k-1 2 k~l 

\ -k -k k-1 2 ) 



5[(/c-2)(/c-l)/c(/c + l) + l] (36) 
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is odd, G5 acts transitively on the space of quadruples {a,b,c,d) (mod 2''^), and 
the subgroup that preserves a given quadruple (and hence the first N rows of each 
horoball in the tiling), is an index 2^^ subgroup of G5, and hence is a cofinite 
subgroup of PSL(2, R). The remainder of the lemma, and the theorem, follow as in 
triangular case. | 

Next we consider the question of conjugacy for these systems. Again, we just 
modify the argument that worked for triangle tilings. 

Theorem 8. Tk^w and T^^uj' are topologically conjugate if and only if w — w' G 5Z. 
%u and %u' are almost conjugate for any w,w'. 

Proof. The proof is essentially the same as the proof of Theorem 5, in particular 
that of almost conjugacy, which we do not discuss further. li w — w' e 5Z, the 
conjugacy is simply adding {w — w')/5 to each dyadic index. The converse follows 
from the analogues of Lemmas 3-6. The proofs of Lemmas 3, 5, and 6 carry over 
almost word-for-word. Lemma 4 was algebraic, and used specific properties of the 
fixed-sum rule for triangles. In its place we have the following two lemmas that are 
specific to pentagonal horoball packings and the Hecke group G5. 

Lemma 8. The only elements of Q[\/2] that appear as elements of matrices in G5 
are —1, and 1. 



1 A 
1 



Proof. Viewed as matrices, the pentagonal Hecke group is generated by 

and (^^i ■ The matrix elements are manifestly elements of Z[A], and 

Z[A] n Q[\/2] = Z. We will show that the only integers that actually appear as 
matrix elements are and ±1. 

Rosen [Rose] showed that an element of Z[A] is an element of a matrix in the 
group if and only if it is (up to sign) the denominator of a finite approximant of the 
continued fraction 

where = ±1 and each is a positive integer (except possibly tq, which may be 
zero). The continued fraction expansion of a real number is not unique, but can 
always be expressed in a unique "reduced form" , one of whose requirements is that if 
Tn — 1, then en+iTn+i 7^—1- The denominators Qn of the successive approximants 
to the (possibly infinite) continued fraction satisfy the recursion: 

Qn = rnXQn-l + €nQn-2, (38) 

and we may take Q-i = and Qo = 1. Writing Q„ = anA -|- bn, the recursion 
becomes: 

an = rn{an-i + bn-l) + enOn-2 (39) 
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(40) 



We claim that the coefficients satisfy three properties: a) an > On-i, (3) hn > ^ 
and 7) an > bn-i- These are easily checked for n = 1,2,3. We prove these 
hold for all n by induction. Suppose they hold for n up to k. We have a^+i = 
fk+iiO'k + bk) + ^k+icik-i- If €k+i = 1 this is manifestly at least ak- If rk+i > 1 and 
Cfc+i = -1 then ak+i - ak ^ rk+ibk + {vk+i -1)0^-0^-1, which is non-negative 
since a/. — ak-i > 0. Finally, if rk+i = 1 and ek+i = —1, then > 2, and by 
property 7 we have 

ak+i = ak + bk-ak-i = ak + {rk - l)ak-i +tkbk-2 > ak + ak-i - bk-2 > ak, (41) 

which is the needed induction for a. Next, bk+i = rk+iak + Ck+ibk-i > ak — bk-i > 
0, which is the needed induction for p. For 7 we note 0^-1-1 ^ ak ~\~ bk — ak—i ^ bk, 
which completes the induction. Finally, since the sequence ak is nondecreasing and 
since ai is positive ak is never zero and Qk is never rational for A; > 1. | 

Lemma 9. A conjugacy between Tk^w and Tk^w' must preserve the locations of 
horoballs exactly. 

Proof. We have already shown that the points on the sphere at infinity where the 
horoballs touch are not changed, that their radii change by a bounded amount, and 
that there are only a finite number of possible values for that change in radius. 
This implies that knowing the first digits of all five indices at a pentagon will 
determine which horoballs grow and shrink, and by how much. Note also that the 
deep interiors of all horoballs are identical, so the change in radius is the same for 
all horoballs, modulo log(2). 

Now consider a tiling whose associated horoball packing is as follows: One 
horoball is the set {x + iy\y > 1}, and the others are its images under the Hecke 



X-axis at Q;/7 with Euclidean diameter 1/7^. 

Note that the dyadic indices of the tiling modulo 2^ are unchanged by the 
transformation z ^ z -\- 2^+^ A, so that (j) of this tiling corresponds to a packing 
that is invariant under z 2+2^+^A. However, if the packing contains the horoball 
{x + iy \y > v}, then it is invariant only under addition of multiples of vX. Thus v 
must divide 2^+^, and in particular must be rational. 

Now consider a horoball that meets the horoball at infinity in the new pack- 
ing (but did not meet the horoball at infinity in the original packing). Since its 
(hyperbolic) radius has changed by — log(j^) (mod log(2)), and since u is rational, 
its Euclidean radius must have been a power of 2 times the square of a rational to 
begin with. However, this implies that there is an element of the Hecke group with 
7 of the form of a product of a rational and a power of ^/2, and not equal to or 
±1. By Lemma 8 no such element exists. | 

Next we consider some differences between optimization problems with other 
variations on our basic tile. So far we have considered packings of by horoballs 
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which meet either in "triangles" (the densest packing of horoballs) or in "pen- 
tagons" , and modified our basic tile to have 2 prongs, each of which is either a third 
of a triangle (Figure 7) or a fifth of a pentagon (Figure fO). One can easily allow 
horoball packings defined by other regular n-gons, and also consider tiles to have 
more prongs, one for each neighboring n-gon; for n = 3 and m = 3 see Figure 11. 

Let T(n, m) denote the space of all tilings by m-pronged tiles, each prong con- 
gruent to one of the n isosceles triangles dividing an n-gon as defined above. Let X 
be a closed, PSL(2, M)-invariant subset of T{n, m) and let y be a closed, PSL(2, M)- 
invariant subset of T{n',m'). We want to show that, under various assumptions, 
X and Y cannot be topologically conjugate. 

Lemma 10. If cf) : X ^ Y is a, topological conjugacy, (f) maps degenerate tilings 
to degenerate tilings and nondegenerate tilings to nondegenerate tilings. Moreover, 
for nondegenerate tilings the points of tangency of the horoballs at the sphere at 
infinity are not changed by (f). 

Proof. The proof is essentially that of Lemma 3. The fact that we were dealing 
with fixed-sum tilings (with n = n' = 3 and m = m' = 2) was never used. 

Theorem 9. If n = 3 and n' = 5 and X C T(n, m) contains nondegenerate tilings 
then X and Y C T{n',m') are not topologically conjugate. 

Proof. For a nondegenerate tiling in T(3, m), the "cusp point set", the set of points 
at infinity that meet horoballs, is conjugate, by some fixed element of PSL(2, M) 
to Q U {oo}. The cusp point set for the corresponding (nondegenerate) tiling in 
T(5, m') is conjugate to Q[A] U oo, where A is the golden mean. This contradicts 
Lemma 10. | 

Theorem 10. U m ^ m' then X C T(n, m) and Y C T(n', m') are not topologi- 
cally conjugate. 

Proof: Assume without loss of generality that m < m' . X necessarily contains all 
the degenerate tilings, and in particular contains a tiling invariant under the map 
z — > mz. By Lemma 10, any conjugacy would have to take this to a degenerate tiling 
in T{n',m') that is also invariant under z mz. However, under this symmetry, 
points along the y-axis are only moved a distance ln(m), while any symmetry of a 
degenerate tiling in T{n' , m') must move points at least a distance In(m'). | 

Note that Theorem 9 requires that X contain nondegenerate tilings; if X con- 
sists only of degenerate tilings, then n is irrelevant. Theorem 10, however, only 
depends on the existence of degenerate tilings. These arise automatically since they 
are in the orbit closure of every nondegenerate tiling. Generalizing Theorem 9 would 
require knowing more than we do about cusp point sets for Hecke groups. 

V. Isomorphism and Uniqueness in Problems of Optimally Dense 
Packings 

Optimally dense packings, especially tilings of and E^, have been impor- 
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tant for many years in classifying certain geometric properties of patterns; we refer 
here to the classification through symmetry groups, the so-called crystallographic 
symmetries. 

At heart the formalism consists of treating the structures of interest as subsets 
of S", with the action on them of Q'^ - that is, one uses the structure of dynamical 
systems. Consider for instance two ball packings in E^, P^'^^, the face centered 
cubic, and P'^^p^ the hexagonal close packed. When we choose to distinguish pf'^^ 
from P^'^P on symmetry grounds what we are saying is that the subgroup of (the 
connected Euclidean group), which acts trivially on every element of one orbit, is 
different from the symmetry group of the elements of the other orbit. This implies 
the systems are not conjugate: there is no bijection, between the two orbits under 
Q^, which intertwines the action of Q^. So in this simple situation we see that 
conjugacy can detect differences in symmetry. 

The case of the kite & dart tilings (Figure 3) is instructive. As we noted in sec- 
tion III, it is natural to want to think of all these tilings as equivalent. This is true 
even though the tilings can actually have different symmetry groups; for instance 
there are two noncongruent kite & dart tilings with a point of 5-fold rotational 
symmetry, which the other kite & dart tilings do not have [Card] . Furthermore, the 
symmetry of the two special tilings actually do not play an essential role, for two 
reasons. First, the 5-fold rotational symmetry appears in regions of arbitrarily large 
size in every tiling, and this could replace the exact symmetry of the special tilings. 
Furthermore the Penrose tilings have a statistical form of 10-fold rotational sym- 
metry, which is expressed by the 10- fold rotational symmetry of all the translation 
invariant measures on the space of Penrose tilings [Radi] . We also note an analogy 
between the different "symmetry" of elements of tilings with different fixed-sums, 
as evidenced by nonconjugacy, and the different (5-fold rotational) symmetry that 
appears among kite & dart tilings. 

So we are led to relax the strict form of equivalence whereby two packings are 
equivalent if they are in the same orbit under Q'^. From the example of the kite 
& dart tilings one might be led to replace this by having optimally dense packings 
equivalent if they are generic for the same measure. But from the various examples 
of section IV we will go one step further. 

Definition 4. We say that two optimally dense packings p' are "weakly equiva- 
lent" if the optimal measures /Up, for which they are generic have the following 
property: the set M{p), of all optimal measures the support of which intersects 
the support of Hp, coincides with the set M{p'). An optimal density problem will 
be said to have a "unique solution" if there is only one weak equivalence class of 
optimally dense packings. 

As we saw in section III, there are simple examples of optimally dense packing 
problems, in particular that for disks of fixed radius in E^, for which the solution is 
unique in the sense of consisting of a single orbit of or, put another way, in the 
sense that the quotient Pq/Q^ consists of a single point. Because of the aperiodicity 
of the kite & dart tilings, and the modified binary tilings, we have been led to divide 
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Pg by a cruder equivalence relation. 

This paper is a continuation of a long tradition of classifying a pattern through 
the dynamical system associated with it by the action of the isometry group of the 
ambient space of the pattern. A common step taken when following the dynamics 
approach is to settle on a form of conjugacy, typically either measurable conjugacy 
or topological conjugacy - and in effect declare two optimization (or tiling) problems 
equivalent if their dynamical systems are conjugate in the chosen sense. Prominent 
in this vein is the analysis by Connes, Putnam, Kellendonk et al noted above, in 
which invariants of aperiodic tilings are sought through operator algebras associated 
with their dynamical systems. 

Definition 5. We declare two optimal density problems, associated with finite sets 
B and B' of bodies in some fixed to be "equivalent" if there is a topological con- 
jugacy between their dynamical systems, Pg,^") and (Pg,^"'), where Pg denotes 
the closure in Pg of Pg. 

In terms of this notion of equivalence the proofs of Theorems 9 and 10 show 
how the optimization problems, for different variations of our modified binary tile, 
can be distinguished by geometric features. 

The operator algebra approach noted above is a powerful way to obtain the de- 
sired invariants for topological conjugacy. Associated to a dynamical system (X, G) 
is the crossed-product C*-algebra C{X) x^G, where C{X) is the C*-algebra of 
continuous complex- valued functions on X and a is the action of G on C{X). (The 
crossed-product is the completion of the algebraic tensor product in a certain norm; 
for this and other terms in operator algebras we refer to [Blac].) The X-theoretic 
invariants of the crossed-product algebra are topological conjugacy invariants for 
the dynamical system. 

A common way to compute invariants is to associate an AF algebra with the 
dynamical system, "large" in an algebra Morita equivalent to the crossed-product 
of interest. (Morita equivalence preserves the K-theoretic invariants.) We do not 
see how to do that here. Alternatively one could try to compute Kq using tools such 
as the Pimsner-Voiculescu 6-term exact sequence; this has been a practical route at 
least when the dynamical group is M or M^, but this seems to be harder for groups 
such as PSL(2,M). 

In short, the operator algebra methods used to produce dynamical invariants for 
aperiodic systems in Euclidean space seem to need extension for this more general 
aperiodic setting. 

VI. Summary 

We have considered in a unified fashion, for Euclidean and hyperbolic spaces, 
the class of problems of determining the optimally dense packings by congruent 
copies of bodies from some fixed finite set B. 

From a qualitative point of view the main solved examples in the literature are 
as follows. For balls of fixed radius in or fixed tight radius in or there is a 
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unique periodic solution [BoRl] . For balls of fixed radius in the problem is widely- 
believed degenerate, with nonunique solution, but including periodic special cases. 
For the general problem for balls of fixed radius the problem is better understood 
in than in E"^, at least in the qualitative sense we are emphasizing; for tight 
radius the problem has the obvious unique periodic solution, and for most radii (all 
but countably many) the problem in is aperiodic, that is, none of the optimally 
dense packings have symmetry group with fundamental domain of finite volume 
[BoRl]. (For non-tight radii none of the optimal packings have been published.) 
For optimal density problems of polyhedra there is again the phenomenon of aperi- 
odicity, exemplified by the kite & dart tilings in E^. Aperiodicity was discovered in 
the Euclidean context by Berger in 1966 [Berg], and has been of growing interest, 
from many perspectives (though mainly still in Euclidean spaces) ever since [Radi] . 

We emphasize how short the above list is, how few optimum density prob- 
lems have been solved. As part of this small but distinguished class, the aperiodic 
ones must play a significant role in any general understanding of optimum density 
packings. 

To prove existence of optimally dense packings it has been found useful to work 
in an ergodic theoretic formalism; in a sense, the solutions of such problems are 
naturally organized through invariant probability measures on a topological space 
of packings. Much of this paper is an analysis within that formalism of certain 
new examples: the densest packings of the hyperbolic plane by congruent copies 
of a certain polygon, the modified binary f, and others in a 2-parameter family 
of variations. These are the first aperiodic optimal density problems in hyperbolic 
space for which explicit solutions have been determined, and they exhibit features 
unknown from examples set in Euclidean spaces - specifically, the small, compact, 
invariant subset T^eg appearing in the supports of all invariant measures on the 
solution space. 

Optimal density problems have a long tradition but have never been treated as 
have other classes of optimization problems - for instance by analysis of conditions 
for existence and uniqueness of solutions. In this paper we consider the question of 
uniqueness, guided by symmetry properties of solutions. We analyze such problems 
as dynamical systems, with the group of isometrics of the ambient space (E"^ or 
W^) acting on the (compact, metrizable) space of all possible packings. We are led 
to classify such systems up to topological conjugacy, and use geometric features as 
invariants. 

More specifically. Theorems 1 and 2 led us to partition the class of optimal 
density problems with unique solution into two classes: the periodic and the ape- 
riodic. The study of periodic optimization problems in E"', for which the solutions 
have cofinite symmetry group, led, many years ago, to classification of the discrete 
subgroups of the isometry group of Euclidean space in low dimensions. Aperiodic 
tilings have led to related work; among problems set in E"^, the study of the quaqua- 
versal aperiodic tilings [CoRa] led to new results on classification of certain (dense) 
subgroups of SO{3) [RaSl], [RaS2], [CoRS]. Similarly, the classification results in 
Section IV of our optimization problems in naturally led to noncofinite sub- 
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groups of PSL(2R), such as the symmetry groups of fixed-sum tihngs (Theorem 6), 
as well as questions about Hecke groups. 

In summary, classification of aperiodic optimization problems amounts to a 
study of the "symmetries" of the packing solutions in senses related to, but dif- 
ferent from, the manner appropriate for the well studied periodic structures. The 
mathematics that is generated by such analysis, a mixture of dynamics, operator 
algebras and Lie groups, is perhaps the main significance of the study of optimal 
density problems. 

Acknowledgments. We are grateful for useful discussions with Alan Reid, in 
particular for pointing us to reference [Rose]. 
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Figure 1. The densest packing of equal disks in the Euchdean plane 
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Figure 2. The Penrose kite & dart tiles of the Euchdean plane 
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Figure 3. A Penrose kite & dart tiling of the Euclidean plane 
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Figure 4. The simple binary tile in the upper half plane model of the hyperbolic plane 



30 



1 1 I I I I f I 



Figure 5. The binary tiling of the upper half plane 
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Figure 6. A packing of the upper half plane by horoballs 
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Figure 7. The modified binary tile in the upper half pL 
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Figure 8. A packing of the Poincare disk by horoballs 
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Figure 9. Another packing of the Poincare disk by horoballs 
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Figure 10. A "pentagonal" binary tile 
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Figure 11. A triangular tile with 3 prongs 
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